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We use the Hirsch-Fye quantum Monte Carlo method to study the single magnetic impurity 
problem in a two-dimensional electron gas with Rashba spin-orbit coupling. We calculate the spin 
susceptibility for various values of spin-orbit coupling, Hubbard interaction, and chemical potential. 
The Kondo temperatures for different parameters are estimated by fitting the universal curves of 
spin susceptibility. We find that the Kondo temperature is almost a linear function of Rashba spin- 
orbit energy when the chemical potential is close to the edge of the conduction band. When the 
chemical potential is far away from the band edge, the Kondo temperature is independent of the 
spin-orbit coupling. These results demonstrate that, for single impurity problem in this system, the 
most important reason to change the Kondo temperature is the divergence of density of states near 
the band edge, and the divergence is induced by the Rashba spin-orbit coupling. 
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Spin-orbit coupling , the interaction between a quan¬ 
tum particle’s spin and its motion, has attracted much 
attention in condensed matter physics. The spin-orbit 
coupling is crucial for quantum spin Hall effect ’ and 
topological insulators 1-11 . Magnetic doping in these sys¬ 
tems is very important. It is reported that magnetic 
doping is an efficient method to tune and to control the 
transport properties of these systems. Magnetic doping 
is a feasible method to break the time reversal symmetry 
on the surface of topological insulators 7-12 , so that one 
may observe very interesting physical phenomena, i.e., 
the topological magnetoelectric effect 13 , the half-integer 
quantum Hall effect 14 , the image magnetic monopole in¬ 
duced by an electric charge near the surface of topolog¬ 
ical insulator ’ , the topologically quantized magneto¬ 

optical Kerr and Faraday 17,18 rotation in units of the 
fine structure constant 11 ’ 20 , the repulsive Casimir inter¬ 
action between topological insulators with opposite topo¬ 
logical magnetoelectric polarizabilities 21-23 , etc. It is also 
demonstrated experimentally and theoretically that mag¬ 
netic doping is a straightforward approach to generate 
the quantum anomalous Hall effect in topological insula¬ 
tor thin films. 

However, as far as we know, the magnetic impurity 
problems, i.e., the Kondo problems 24 , in spin-orbit cou¬ 
pled systems are far from clear up to now. A clear un¬ 
derstanding of the Kondo problem is essential to study 
the low temperature transport properties and the mag¬ 
netic ordering problems in these systems. In this paper, 
we study the most fundamental quantity of the Kondo 
problems, the Kondo temperature, in two dimensional 
electron gas with Rashba spin-orbit coupling. There 
are some previous works have studied the Kondo tem¬ 
perature of spin-orbit coupled systems 2 - . A unitary 
transformation is applied to demonstrate that the spin- 
orbit coupling can be absorbed into the effective hopping 


amplitudes and it has no effect on the Kondo tempera¬ 
ture. Mean-field calculations based on variational wave 
functions show' that the Kondo temperature is expo¬ 
nentially altered by the spin-orbit coupling, because the 
density of state on the Fermi surface is changed by spin- 
orbit coupling. Perturbation renormalization group anal¬ 
ysis based on the Schrieffer-Wolff transformation shows 
that the Kondo temperature is exponentially enhanced 
by the Dzyaloshinskii-Moriya (DM) interaction between 
the local moment and the itinerant electrons. The numer¬ 
ical renormalization group (NRG) calculations show 28 
that the Kondo temperature is almost a linear function 
of Rashba spin-orbit coupling energy. A dependence of 
the Kondo temperature on the spin-orbit coupling is also 
expected for quantum dots 2 ”, where due to additional 
parameters like left-right asymmetries in the tunneling 
couplings or magnetic fields more complicated behaviour 
including suppression may appear. These qualitatively 
different results show that the variation of Kondo tem¬ 
perature as a function of spin-orbit coupling in this sys¬ 
tem is still an open question. 

In this work, we use the Hirsch-Fye quantum 
Monte Carlo 11 (HFQMC) simulation to study the prob¬ 
lem. HFQMC is a numerically exact method. It 
has been widely used to study the magnetic impu¬ 
rity problems in normal metals, i.e., the local mo¬ 
ment of magnetic impurity 30,31 , the Kondo effect 2 ’ 2 , 
the Ruderman-Kittel-Kasuya-Yosida interaction between 
magnetic impurities 33 , etc. Recently the HFQMC tech¬ 
nique has been applied to study the local moment 
formation of Anderson impurities in dilute magnetic 
semiconductors and graphene ’ 5 , and the interaction 
between magnetic impurities in the bulk of topological 
insulators 2 ' 1 . 

Let us present the model Hamiltonian here, the total 
Hamiltonian has three different terms: the bulk state 


2 


of two-dimensional electron gas with Rashba spin-orbit 
coupling 1 ’''’ 28 , the magnetic impurity with strong on¬ 
site Coulomb interaction, and the hybridization between 
the electron gas and the impurity states, 


Ti. = T~Lh + kt d + 'Hhyb, 


(1) 


T~Lb = Sfc, s =t,4. <S 2m M d" a(k x a y — k y a x ) i/j ks ,(2) 
kid = 4( e d - p)d s + Ud\d^d[di, (3) 

"Hhyb = ^/fc,s=T 4 {Vk1pk,sd s + V k rfg'i/’fc.s) , (4) 


where ip k and 'tp kyS are creation and annihilation oper¬ 
ators of two-dimensional electron gas with momentum k 
(the Plank constant h has been set to be 1) and spin s. 
m is the effective mass of two-dimensional electron gas. 
p is the chemical potential, a is the strength of Rashba 
spin-orbit coupling. cr x , y , z are the spin Pauli matrices. 

(dj) and d-j- (djJ are creation and annihilation opera¬ 
tors of spin-up (spin-down) impurity states, respectively. 
e d is the energy level of the impurity state. U is the 
on-site Coulomb interaction. V k is the hybridization be¬ 
tween bulk electrons and impurity states. Without lose 
of generality, we set the hybridization to be short range 
coupled, so that the hybridization Vk is independent of 
momentum k. 

Here we make a note on the units and parameters 
used in this paper. In the continuum limit, the summa¬ 
tion over momentum is approximated by an integration, 
Efc = r ( 2 J) 2 i w here N 0 is the number of sites, 

kx is determined by the band width D = k^/2m, the 
total density of states (DOS) is p = No/2D , and con- 
verntionlly we define the hybridization T = npV 2 . We 
set r = 0.057T and U = 1.0 in numerical calculation so 
that our results can be comparable with those given in 
numerical renormalization group (NRG) calculations. 

The bulk states can be integrated out to get the par¬ 
tition function of impurity states, 

Z = J Vd\Vd s e~ s (5) 

5 = dr {dt [£ + (e d — p) + E d ] d s + C/d{d t d}d ; }(6) 

where /3 = 1/ksT is the inverse temperature, ks is the 
Boltzmann constant. And is the self-energy of d- 
electrons, which is given in imaginary-frequency repre¬ 
sentation by, 




FIG. 1. (a) The product of temperature T and spin suscepti¬ 
bility x as a function of temperature for different values of U. 
The spin-orbit coupling is chosen to be E a = ma 2 /2 = 0.5. 
Chemical potential p = 0.2, and impurity energy td — p = 
—U/2. (b) The universal curves of T\ for rescaled temper¬ 
atures (see the context for more details about rescaled tem¬ 
peratures) . 


quantities, i.e., the local charge, the local moment, the 
spin susceptibility, etc., from the imaginary-time Green’s 
function. The dynamical quantities, i.e., the spectral 
function, can be obtained in HFQMC by the maximal 
entropy method ’ . The spin susceptibility x is an i m_ 
portant quantity to describe the Kondo problem, 

X= [ dT(S z (T)S z ( 0)). 

Jo 


Ed(iwn) = — log ( - + ) + 2Carctanh( + 7 

7T |_ \ 47 T J 

( 7 ) 

where C, = a/^/a 2 + 2 (iui n + p ), a = y/ma, and 7 = 
0.5772 is the Euler-Mascheroni constant, from dimen¬ 
sional regularization(Make a note here: why I need the 
dimensional regularization). The DOS of bulk electrons 
is proportional to the imaginary part of T, d {oj + iS). 

In HFQMC simulation, the Hubbard interaction is sim¬ 
ulated by the statistical average over auxiliary Ising spin 
configurations. One can calculate many important static 


where S z = d|d-|- — dj)dj,. At very high temperature (i.e., 
T > 10 in Fig. 1(a)), the total system behaves like 
free electron gas, which satisfies the Curie-Wiess law and 
T\ = 1/2 ;! . At the intermediate region (0.1 < T < 1 
in Fig. 1(a)), the impurity behaves as a local moment 
state. Our numerical simulation show that the square 
of the local moments m 2 for the parameters given in 
Fig. 1 are enhanced in this intermediate region. At very 
low temperature, the local moment state is screened by 
the itinerant electrons, and the product of temperature 
and spin susceptibility satisfies the following universal 
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FIG. 2. The product T\ for different impurity energies and 
different spin-orbit couplings, (a), (c) and (e) are the original 
data for different impurity energies Ed — p = —0.3, Ed — p = 
—0.5 and £d~ P = —0.7 respectively, (b), (d) and (f) are the 
corresponding rescaled curves of (a), (c) and (e). Hubbard 
interaction U = 1.0 and chemical potential p = 0.2 are chosen. 


relationship , 

$(4Tx(T)-l) = ln(T/T ir ), (8) 

where $ is a universal function and Tk is the Kondo tem¬ 
perature. We can use this relationship to fit the variation 
of the Kondo temperature as a function of different pa¬ 
rameters, i.e. Hubbard interaction, spin-orbit coupling, 
chemical potential, hybridization, etc. 

Before a more detailed discussion about the effect of 
spin-orbit coupling, we show the universal curves of spin 
susceptibility for different Hubbard interactions, shown 
in Fig. 1. One can see that the curves in Fig. 1(a) 
are almost equal-spacing under the logarithmic coordi¬ 
nate of temperature, which demonstrates that the Kondo 
temperature is exponentially altered by the Hubbard in¬ 
teraction. We can use the following formula to fit the 
Kondo temperature, 

T k (U) =T ( ° ] e~ aU . (9) 

Fig. 1(b) is the numerical fitting of the universal relation¬ 
ship Eq.(8) about Fig. 1(a), where the horizontal coordi¬ 
nate is fixed and the longitudinal coordinate is rescaled 



FIG. 3. Kondo temperature as a function of spin-orbit cou¬ 
pling for different impurity energies. Hubbard interaction 
U = 1.0 and chemical potential p = 0.2. 


for different Hubbard interactions. We find that the fit¬ 
ting result a = 1.11572. This value is different from the 
one estimated from the ’’flat band” results 

T k =D^Jex P (-l/pJ) (10) 

where pJ = ST/nU and D is the effective bandwidth. Eq. 
(10) gives that a = 2.5 for the hybridization T = 0.057T 
used here. This result demonstrates that the density of 
states near the band edge is important to alter the Kondo 
temperature. 

Now we study the Kondo temperature for different 
spin-orbit couplings. There are three important ap¬ 
proaches which may change the Kondo temperature: (1) 
The DM interaction between the local moment of the 
impurity and the itinerant electrons, which may expon- 
tionally enhance the Kondo temperature; (2) The spin- 
orbit coupling can spit the two-fold spin degeneracy of 
Fermi surface, so that the DOS on the two Fermi surface 
may change the Kondo temperature; And (3), although 
the total DOS on the Fermi surface is not altered by the 
spin-orbit coupling (see Fig. (4) for more details), the 
divergence of DOS near the band edge may change the 
Kondo temperature. We find that our results support 
case (3). Now we show the numerical results. 

Fig. 2 gives the product of temperature and spin sus¬ 
ceptibility as a function of temperature for different spin- 
orbit couplings and different impurity energy levels. Figs. 
2 (a), 2 (c) and 2 (f) are original data for impurity energies 
ed — /i = —0.3, Ed — p = —0.5 and — p = —0.7 re¬ 
spectively. Figs. 2(b), 2(d) and 2(f) are the rescaled 
curves. The numerical fitting of the Kondo tempera¬ 
tures are given in Fig. 3. We find that our QMC re¬ 
sults are qualitatively consistent with those from NRG 
calculations. The Kondo temperatures are almost lin¬ 
ear functions of the Rashba energy E a = ma 2 /2. How¬ 
ever, there are some differences, ( 1 ) when the impurity 
energy e — p = —0.5 = —17/2, the NRG calculations 
show that the Kondo temperature is enhanced by the 
spin-orbit coupling, our QMC simulations show that it is 
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FIG. 4. (a) The product T\ for different spin-orbit couplings 
when the chemical potential p = 3.0, Hubbard interaction 
U = 1.0 and impurity energy ea — p = —U/2. (b) The DOS 
of the bulk electrons for different spin-orbit couplings. The 
solid lines are the total DOS, dash lines and dotted-dash lines 
are the DOS of the two spitted bands. 


suppressed; (2) the effect of spin-orbit coupling in QMC 
simulation is much smaller than those given in NRG. 
These differences demonstrate that different identifica¬ 
tions of the Kondo temperature may have some quanti¬ 
tative differences. More importantly, one can find that, 
when the temperature is in the intermediate region ( i.e. 
0.1 < T < 1), all of the three cases (e^ — p = —0.3, —0.5, 
— .07) have the tendency that the Kondo temperature is 
enhanced by the spin-orbit coupling. This demonstrates 
that maybe the perturbative renormalisation analysis is 
valid in this parameter region. However, when the tem¬ 
perature come into the Kondo region, i.e. T < 0.1, the 
perturbative renormalisation analysis breaks down, and 


there may exist some crossover from the enhancement to 
the depression, i.e. T ss 0.1 in Fig. 2(e). 

These analysis show that the effect of DM interaction 
can be neglected at very low temperature. There are 
two other important effects which can change the Kondo 
temperature: the difference of the two Fermi surfaces and 
the divergence of the DOS near the band edge. Now we 
discuss the results when the chemical potential is far away 
from the band edge. The numerical results are shown 
in Fig. 4, we calculate the product of temperature and 
the spin susceptibility as a function of temperature for 
various values of spin-orbit coupling from 0.1 to 3.0. In 
Fig. 4(b), we can see that the difference between the DOS 
of the two bands is dramatically enhanced. However, the 
universal curves in Fig. 4(a) are slightly changed and 
the Kondo temperature is almost unchanged by the spin- 
orbit coupling. 

So we conclude that the exponential enhancement of 
the Kondo temperature from perturbative renormalisa¬ 
tion group analysis of the DM interaction between the 
impurity state and the bulk electrons breaks down in the 
low temperature Kondo region. Kondo temperature is al¬ 
most a linear function of Rashba spin-orbit energy when 
the chemical potential is close to the edge of the conduc¬ 
tion band, and when the chemical potential is far away 
from the band edge, the Kondo temperature is indepen¬ 
dent of the spin-orbit coupling. And the divergence of the 
DOS near the band edge is the most important factor to 
alter the Kondo temperature. 
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